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Abstract 

A lepton flavor violating process \i — > is investigated in the supersymmetric extra 
U(l) models, which often appear as the low energy effective models of superstring and 
can potentially solve the /z-problem. The branching ratio of this process is calculated. It 
is numerically estimated and compared with that of the MSSM. In this study we take 
account of an abelian gaugino kinetic term mixing and discuss its influence on this process. 
The possibility to find the extra gauge structure through this process is discussed. 
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The standard model(SM) has shown its incredible accuracy to describe the electroweak 
interaction through the precise measurement at LEP. Nevertheless, physics beyond the 
SM is eagerly explored because of its unsatisfactory feature for the explanation of the 
origin of weak scale and its stability. Supersymmetrization of the SM is now considered 
as the most promising extension to solve this problem]]]]. However, even in this minimal 
supersymmetric standard model(MSSM) there still remains a theoretically unsatisfactory 
feature which is known as the /i-problemf2|]. To cause an appropriate radiative symmetry 

I/O 

breaking at the weak scale, we need a Higgs mixing term /zifii^, where /i ~ 0(G F ±/Z ) and 
Gf is a Fermi constant. However, there is no reason why /i should be such a scale because 
it is usually considered to be irrelevant to the supersymmetry breaking. A solution for 
this problem is to consider /i as a dynamical variable ||. The introduction of a singlet 
field 5* with a Yukawa type coupling XSH1H2 can realize this scenario in the simplest 
wayj|]. That is, if 5* gets a vacuum expectation value(VEV) of order 1 TeV as a result 
of radiative corrections to the soft supersymmetry breaking parameters ||, /i ~ 0(G F 1 ^ 2 ) 
will be realized dynamically through the relation /i = A (S). 

The extra U(l) models are the typical extensions of gauge structure of the SM. It is 
very interesting to note that many extra U(l) models have the above mentioned feature 
inevitably || [7]]. Low energy models derived from superstring often have accompanied 
with extra U(l) factors in their gauge structure ||. It seems to be natural that these 
aspects motivate us to investigate extra U(l) models and try to look for a clue of such a 
gauge structure. Recent precise measurements at LEP and also the Tevatron experiments 
show us that the lower bound for its gauge boson mass is rather large and then it may 
not be so easy to find it directly ||. Even in that case if nature is supersymmetric and the 
gauge bosons have their superpartners, there may be other possibilities to investigate the 
gauge structure through examining the processes to which their superpartners contribute. 

In this letter we study the lepton flavor violating /i — > process. The gauginos of 
extra U(l)s can affect this process. Our purpose here is to estimate their effect and discuss 
the possibility to find the extra gauge structure through this process. Its comparison 
with the results in the MSSM will also be useful for the future experimental analysis. We 
consider the minimal models which have only one extra U(l)x and a singlet Higgs S with 
a U(l)x charge besides the MSSM contents.[] These fields are assumed to remain light 

1 In order to induce the symmetry breaking radiatively, it is necessary to introduce the vector like 
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around the TeV region. The neutralino sector in this model is extended by an extra U(l)x 
gaugino and a fermionic partner of the singlet Higgs 5" in addition to the ingredients of 
the MSSM. 

Before proceeding to the detailed study we should note an additional feature of the 
neutralino sector of these extra U(l) models. It has been well known that in principle 
there can be kinetic term mixings among abelian gauge fields because these field strengthes 
are gauge invariant .0 Super symmetrization of the models introduces kinetic term mixings 
among abelian gauginos. In the analysis of multi U(l) models we generally need to take 
account of these effects. Thus at first we briefly summarize the mixing effects in the 
gaugino sector for the usage in the later study. 

In supersymmetric models gauge fields are extended to vector superfields 

v wz {x, e, o) = -ea^ev^ + ieeex - mex + -eeeeD, (i) 

where we used the Wess-Zumino gauge. A gauge field strength is included in the chiral 
superfield constructed from Vwz in the well known procedure, 



4zA Q - A6 a D + U9Pa,<>%{d»V v - PV>) - Wa^xK (2) 



W a {x,0) = (DD)D a V wz 
In terms of these superfields the supersymmetric Lagrangian can be written as 

c = h ( waw °)F + ( $t eM2g°QVwz)®) D , (3) 

where $ = (<f),ip,F) is the chiral superfield representing matter fields. This Lagrangian 
is easily extended to multi U(l) models. In the models with two U(l) factor groups, the 
supersymmetric gauge invariant kinetic terms are most generally written by using chiral 



superfields W% and W° for U(l) a x U(l) b as[|TJ, 



— (w aa WZ) + — (w ba W<>) +*-^(w aa W h ^ (4) 
32 V a ) f 32 V a J f 16 V a J f v ' 



extra color triplets (g,g) which have the coupling to the singlet S as AcS'ggl^. But in the present study 

they play no role and then we will ignore them. 

2 In fact there are some works in which it is discussed in what case kinetic term mixings can occur |fl0|, 

in. 
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where we introduced the mixing terms. These can be canonically diagonalized by per- 
forming the transformation, 



f W a 



( 



( 



w a 
w b 



(5) 



1 — tanx 
^ 1/cosx J \ 

This transformation affects not only the gauge vector fields but also the sector of gauginos 
A 0> 6 and auxiliary fields D a ^. The modification due to this transformation in the gaugino 
sector can be summarized as 

9°aQa\ a + 9° b Q b X b = 9aQaX a + (9abQa + 9bQb) X\ (6) 

where X a ,b are canonically normalized gauginos. Q a and Qb stand for the charges of 
matter fields for U(l) a and U(l)b- Gauge coupling constants g a , g a b and <?& are related to 
the original ones as, 

9° 

9a = 9° a , 9ab = 9a^X, 9b = ■ (7) 

cosx 

These low energy values are determined by using the renormalization group equations. 
However, in the present study we will treat them as parameters. 

For the study of the /i — > ej process in the supersymmetric models, it is necessary to 
clarify both of the neutralino and chargino sector. The relevant part of the superpotential 
and soft supersymmetry breaking terms are 

W = XSH 1 H 2 + ---, 

i 

+ ^(M W X W X W + MyXyXy + M X X X X X + MyxXyXx + h.c.) + ■ • • , (8) 

where fa represents the scalar components contained in the models. M^, My and 
M x are soft supersymmetry breaking gaugino massesQ for SU(2)l, U(l)y and U{l) x . 
These parameters and a Yukawa coupling A are assumed to be real and it should not 
be confused with the gaugino fields A a . Using the canonically normalized basis, we 
can write down the relevant quantities in the neutralino sector modified by the kinetic 
3 We introduce the effect caused from the abelian gaugino mass mixing as Myx, which may exist at 
the Planck and may also be yielded through the loop effects. We need to estimate its low energy value by 
using the renormalization group equations. In the later numerical study we put Myx = 0, for simplicity. 
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term mixing. They are the neutralino mass matrix and the vertex factors of gaugino- 
fermion-sfermion interactions. If we take the canonically normalized gaugino basis as 
M T = (— i\w 3 , —i^Y, —i^x, Hi,H 2 , S) and define the mass terms as 

the 6x6 neutralino mass matrix M. can be expressed as 
/ 



(9) 



V 



M w 








TJlzCw cos (3 


—m z cw sin /? 
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Ci 


—m z sw cos (3 


m 2 % sin/3 








Ci 
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c 3 
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mzcw cos (3 


—mzsw cos (3 
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Xv sin (3 


—m z c w sin (3 


m z s w sin (3 


c A 


Xu 





Xv cos (3 








c 5 


Xv sin (3 


At> cos f3 






/ 



Matrix elements C\ ~ C5 are components which are affected by the kinetic term mixing. 
They are represented as 



d = -My tan x + 



My , 2 , M x 2M yx sin X 
, 6 2 = My tan x + T7"x 



cosx 



cos^ x 



cos 2 x 



1 / , . #x<2i\ fl _ if. , #x<2 2 \ . fl 

-7^ 9y tan x H vcos/3, C 4 = -7= — #y tan x H usm/3, 

V2 v cos % y V2 V cos X J 



n 1 9xQs 



(10) 



y/2 COSX 

where Q±, Q2 and are the extra U(l)x charges of Higgs chiral superfields Hi, H 2 and 
S. 

Neutralino mass eigenstates Xi (i = 1 ~ 6) are related to Hj through the mixing matrix 

U as 

x? = EW (ii) 

The change induced by the kinetic term mixing in the gaugino interactions can be confined 
into the extra U(l)x gaugino sector and by using eq.(6) new interaction terms can be 
expressed as, 



r (-g Y Ytan X + ?^-) X x 4> 
V cosx / 



1 

+H* ( - 9Y Y tan X + 9j ^A X X H - (-g Y Yt^ X + ) \.n H H 



-g Y Y tan X + 9x ® x j X x ipip 
cosx J 



cosx 



cosx 



(12) 



where ip an d i> represent the quarks/leptons and the squarks/sleptons, respectively. Higgs 
fields (Hi, H 2 , S) are summarized as H and the corresponding Higgsinos Hi, H 2 and 5 
are denoted as H. Taking account of this, gaugino-fermion-sfermion vertecies in the basis 
of mass eigenstates are assigned by the following factors, 

gxQx 



zt(Y,Q 



x, 



1 

71 



gwT 3 Uii + g Y YU 2i + -gyFtanx + 



U: 



cosx 



3/ 



Z*(Y,Q x ) = -j= 



g Y YU 2i + -g Y Y tanx + 



gxQ 



X 



cosx 



U; 



:-Vi 



(13) 



where we used the left handed basis for the chiral superfields. It is also necessary to define 
the chargino mass eigenstates for the following calculation. The chargino mass term is 
given in the matrix form as 

/ 



C 



c 

mass 



— \u \^2mzCw cos (3 

y/2m z c w sin (3 M w 



h.c. 



-i\ 



(14) 



Y 



The mass eigenstates are defined in terms of the weak interaction eigenstates through 
unitary transformations, 



y X2 



Hi 





(15) 



Based on these preparations, we proceed to the estimation of the /i — > process 
in the present models. The flavor changing processes are strongly suppressed through 
the experimental result s|T!2|]. In the supersymmetric models, however, there are generally 
many sources for these processes in the superpartner sector besides the ones of the SM[13- 
16]. Colored superpartners cause the dominant contributions in many hadronic flavor 
changing neutral processes. Thus the contribution from the neutralino sector may not 
be clearly seen through such processes. In order to see the structure of the neutralino 
sector, we need non-hadronic process and /i — > seems to be particurally interesting 
in the relation to our present purpose as far as R-parity violating terms are absent. At 
one-loop level this process can occur because of the existence of nontrivial flavor structure 
of soft supersymmetry breaking terms in the slepton sector. Various studies of /i — ► ej in 
the MSSM framework and some extended models have been done by now[|il| [15 , I7|. We 
extend these analyses to the multi U(l)s case. One-loop diagrams contributing to this 
process in the present models are shown in Fig.l. 
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The effective interaction describing this decay is given as 



(16) 



By carrying out the calculation of diagrams in Fig.l, we can obtain the effective couplings 
Gl and Gr- In the diagonalizing basis of the lepton mass matrix mi, the origin of flavor 
changings is the off-diagonal elements of Kobayashi-Maskawa matrix in the lepton sector 4 
and also the slepton mass matrix. The slepton mass matrix is written asf] 

/ 



ML 



M 2 LR 



M 2 RL M 2 RR 



(17) 



where M\ R = mi(A + Xu tan/?) for charged sleptons. In the neutrino sector right-handed 
sneutrinos are assumed to have the large supersymmetric masses and their relevant part 
to the present analysis is only the part of M\ L . To reduce the number of free parameters 
we make the following assumptions for components of the slepton mass matrix, 



{Ml 



e,u\2 
LL / act 



(m; 



RR) 



M 



(M 



e,u\2 
LLJafS 



(M; 



RR)a/3 



A 2 



As shown in Fig.l, there are two types of diagrams which are distinguished by the place of 
the chirality flip. For our present purpose, it will be enough to concentrate our attention 
on the diagrams with the chirality flip on the internal line ( Figs (a) and (c)).[] Under 
these assumptions the effective couplings Gl and Gr can be summarized as 

(M| ^ 2 t - 2 
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f (MRL)le 


{ M 4 
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c2 



M 2 



A 2 M 2 

—^'tept ( ' ■ 



Ct=fl,T 



m. 



c2 ■ 



Zi(-l,Qe L )Z?(2,Q eR )+ 9wm » m 



u 2 

2mLcos 2 /3 4 \ 



4 We assume the non-zero Majorana neutrino masses induced from the seesaw mechanism in view of 

the solar neutrino problem. 

5 The origin of flavor violating off-diagonal elements of M^ L , M^ R and M\ R are discussed from various 

view points |fl6[ p~7| . In this analysis we donot refer to it and only treat them as parameters. 

6 This treatment may not be bad even in the quantitative view point since the neutralino masses are 

expected to be much larger than charged lepton masses. For the completeness of our formulus, however, 

we will present the contribution to the effective couplings from Figs, (b) and (d) in the appendix. 



7 



9w 



V2m w cos (3 M 4 



(m^(-l, Q e J + m e Zf (2, Q e j) Fi 



, m. 



»,2 



M 2 



r r9 



V2m w cos /3 § (5"" e J( M 2 ) + ^ 



A 2 



M 2 



a=e,r 



m. 



(+)u/(-) 

i 



(19) 



where mf and represent the i-th mass eigenvalues of neutralinos and charginos, re- 
spectively. A 2 e stands for an off-diagonal element between the e- and /^-generation of 
slepton mass matrices as defined by eq.(18). Its allowed range may be estimated at a 



few GeV or less depending on other soft supersymmetry breaking parameters |Lo|. K a p 
is the Kobayashi-Maskawa matrix element in the lepton sector. Kinematical functions 
F\{r) and J(r) appearing from the loop integrals are defined by 

1 



F(r) 
J(r) 



2(1 -r) 4 
1 

2(1 -r) 2 



1 + 4r - 5r + 2r(r + 2) lnr 



-3 + r 



lnr 



Using these results, we can represent the branching ratio of this decay process as 

48tt 2 



(20) 



\Ql\ 2 + \G R ? 



(21) 



In order to compare this result with the MSSM one, it is useful to list up the extra 
parameters added to the ones contained in the MSSM formulus: 

tan/3, (M e LR ) 2 a/3 , M\ A 2 Q/3 , K af3 , M w , My. 

Additional parameters to these are new gaugino masses (Mx, Myx), the kinetic term mix- 
ing parameter sin x, the extra U (1) coupling gx and charges[] (Qi, Q2) and also the /i-term 
relevant parameters (A, (S)(= u)). We can easily check that the neutralino contribution 
to eq.(21) results in the expression given in refs.[fI4, [H| in the case of the photino domi- 



nated neutralino, if we put these additional parameters zero instead of keeping /a(= Xu) 
constant. 

Before choosing a parameter set for the numerical analysis, we should note some 
features of our models. In these models the vacuum expectation value u of the singlet 
Higgs S is relevant to the extra U(l)x gauge boson mass besides determining the //-scale. 
The mixing between the ordinary Z° and the U(l)x boson is severely constrained by 
7 It should be noted that Qi + Q2 + Qs = is satisfied because of the form of superpotential. 
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the precise measurement at LEP and the direct search at Tevatron|§. This constraint 
requires that the mass of the U(l)x boson is large enoughF] and in that case its mass 
eigenvalue is given by[pf 



rn 



,2 _ 1 „2 /^2„.2 , /-|2„.2 /V2 .2 



z 2cos 2 x 



9x(Qi v i + Q2 v 2 + Qs u )- (22) 



The experimental bound on m|, determines the lower bound on u. On the other hand, 
u determines the /i-scale as n — Xu. Thus to keep \x in the suitable range we need to 
put the upper bound on A. For its rough estimation, we take mz> ^ 400 GeV and also 
assume gx = 9y, which is satisfied, for example, in the abelian subgroup of E$ if the full 
components of 27 of Eq contribute to /3-functions. In this case if we require \i < 1 TeV, 
we obtain the upper bound on A as A ^ O.QQs- This bound seems to be reasonable from 
the view point of the analysis of the radiative symmetry breaking |J. 

Another feature which we should note is the dependence of the effective couplings 
Ql and Qr on the neutralino and chargino mass eigenvalues m™ and m?. These de- 
pendence can be factorized as r 1//2 Fi(r) where r = (m"/M) 2 for neutralinos and also 
(r 1 ), r l ' 2 J(r) where r = (m^/M) 2 for charginos. They vary in the range < 
r 1/2 Fi(r) < 0.1, < r -3 / 2 Fi(r -1 ) < 0.25 and < r 1/2 J(r) < 0.44. Each maximum value 
is realized at r ~ 0.27, r ~ 0.025 and r ~ 0.12, respectively. Thus all of these factors can 
be considered as the same order at least except for r ~ 0. 

Taking account of this, we can roughly estimate the condition for the neutralino con- 
tribution dominance by comparing the neutralino contribution to the branching ratio with 
the chargino contribution. The couplings of the neutralinos to leptons come from gauge 
couplings and Yukawa couplings. Because Yukawa couplings are small enough, the dom- 
inant contribution will be yielded by the neutralinos which are dominantly composed of 
the gauginos Xw, Ay and Ax- As seen from eq.(19), it is naively required 



^ue > rn w cos(3 



10 3 (23) 



i. M LR)le ~ rn^gw 

in order to guarantee the similar order contribution from all neutralinos. This requirement 

seems to be difficult to be satisfied. So we focus our study to the case where main 

neutralino contribution comes from the term associated with (M| R ) 2 e in eq.(19). Next 

8 The mixing component of the mass matrix can be small enough for the special value of tan (3 and 
sin x- In such a case this requirement is not necessary to be satisfied. 
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we compare it with the chargino contribution. If we pay attention on the factors in each 
terms of eq.(19) except for the mixing matrix elements, the condision for the neutralino 
contribution becoming larger than the chargino one can be roughly estimated as 

(^-s^^- 20 ^- (24) 

In this estimation we assumed tan/3 ~ 1 and M ~ 100 GeV. If we note that (MJ^) 2 = 
m^A^e and m^tan/3 does not contribute to it, the above condition for (Ml R )^ e corre- 
sponds to A^ > 200-ff^e GeV. Thus if we take K v e ~ 5 x 10~ 4 as the KM-matrix element 
in the lepton sector^] the neutralino contribution assumed above can be expected to be 
dominant under the condition^] 

A Me < 1 GeV, A^ e > 10- 1 GeV. (25) 

This argument suggests that the gaugino components of the neutralino contribution can be 
dominant one in the rather general situation. Moreover, under this condition B([i — > ej) 
takes the value just below the present experimental bound. Thus in such a parameter 
range we may have a chance to get a hint for an additional abelian gauge structure through 
the \x — ► ej process. Our main interest is the effect coming from new ingredients so that 
in the following numerical study we assume the measure of flavor mixing {M^^) 2 so as 
B(fi —>■ &y) to be the same order as the present experimental bound in the MSSM case. 

Now we give our result of the numerical analysis. As the typical values of free param- 
eters, we take 



tan (3 = 1.5, A^ = 0.2 GeV, M = 100 GeV, 
5 

M Y = M X = - tan 2 6 W M W , M YX = 0, A = 0.5, (26) 
3 

where we assumed the unification relation for the gaugino masses. As a typical example of 
the low energy extra U(l), we take the ^-model induced from Eq. Their charge assignment 
for the relevant fields is listed in Table 1. Under this parameter setting, in Fig. 2 B(fi — > ej) 
in this model is drawn as a function of u in the case of Mjy = 80, 180 GeV and 
sinx = 0, 0.3. The horizontal axis should be converted to u by u = 50(u' + 2). Thus the 

9 This small value does not contradict the neutrino oscillation solution for the solar neutrino problem, 
if we assume the existence of a sterile neutrino as ref. |l8| . 

10 These values are very similar to ones given in ref. p5[ . 
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Table 1 The charge assignments of extra U(l)s induced from Eq. These charges are 

normalized as Qi = 20. Only relevant fields to our study are listed from 27 of Eq. 
ie27 

range of u is 100 GeV < u < 2000 GeV. Here it is useful to recall again that mz> is given 
by eq.(22) and also the relation /i = \u. From the recent chargino and neutralino search 
at LEP[|l!|, the allowed region of (/z, M w ) plane should satisfy /z, M w > 100 GeV as far 
as > 0. This corresponds to u' > 2 for A = 0.5. Moreover, taking account of the small 
mixing contraint on Z° and the extra U(l), u > 800 GeV and then u' > 14 should be 
satisfied if we require mz> > 400 GeV,. 

The ratio R of this B(fi — > cy) against the one of the MSSM is also presented in Fig. 3 
in the case of sinx = and 0.3 for each value of Myj/. As easily seen from Figs. 2 and 
3, the influence appearing in B(fi — > ej) due to the extra U(l) gaugino is not so large 
but non- negligible. These figures also show that the gaugino kinetic term mixing has an 
enhancement effect on ^ e^y) in the 77-model. This effect is considered to be mainly 
caused from the change of the effective couplings Z^(Y, Q x ) and ZfiY, Qx)- The smaller 
My/ and u give the larger B(fi — > ej). The effect of sin x 7^ can be easily seen for the 
smaller Mw and u. 

It should be noted that there can remain the deviation of O(10~ 12 ) from the MSSM 
value at u' ~ 14. As u' becomes larger, B(fi — > cy) monotonically decreases but the 
O(10~ 13 ) deviation can still remain. For the smaller u(u' < 14), the suitable condition 
should be satisfied for the consistency with the bound of the mixing between Z° and the 
extra U(l)x, as already remarked in the footnote, if such a condition is satisfied and then 
the small value of u is allowed from the precise measurement at LEP, we may have an 
important clue to study the extra gauge structure because of the large deviation from the 
MSSM prediction. Anyway, these results seem to be interesting. Particularly, the fact 
that B({i —>■ cy) can deviate by O(10" 12 ~" 13 ) from the MSSM value even at the rather 
large u region seems to be encouraging. It may be possible to find some clues of the extra 
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U(l) gauge structure through B(fj, — > ej) if its experimental bound is improved by an 
order from the present value. 

Some comments should be ordered on the various parameter dependences of £>(/i — > 
ej). In particular, A and tan/? dependence seems to be important. 

1. In the present models A can be an independent parameter and it affects B(fi — > cy) 
through the neutralino mass matrix directly besides through the combination with u. 
As a result, B(fi — > ej) shows the substantial A dependence, although no significant A 
dependence can be seen in the ratio R at the large u region. 

2. The tan/3 dependence seems to be very crucial for the absolute value of B{ji — > ej) in 
the whole range of u. The large tan/3 makes B(fi — > e^y) small for all u region and also 
increases the sensitivity to u at the small u region. No significant tan/3 dependence can 
be seen in R at the large u region. 

We may also be interested in the Myx dependence in the case of sin \ 7^ 0. As far as Myx 
is induced by the loop effects, M YX may be roughly estimated as M YX ~ 9y9yxM y /871 2 , 
which is completely negligible. Unless M Y x is produced at M pl as a rather large value, it 
can be safely neglected in this type of study. 

In summary we investigated the fi — > e'j process in the extra U(l) models taking 
account of the gaugino kinetic term mixing. After driving the formulus for the branching 
ratio of fi — > e7 in the general framework, we practiced the numerical study and showed 
that in the ^-model the deviation from the MSSM can be seen at the level of O(10~ 12 ~~ 13 ) 
through this process. The abelian gaugino kinetic term mixing has some effects on this 
process and we may find the suitable clue of the extra gauge structure by investigating this 
process. It will be useful to note that the usage of this process may open the alternative 
window to search the extra U(l) gauge structure for an appropriate parameter region. It 
will be necessary to investigate this process in more general parameter region and other 
extra U(l) models. The combined study with other rare process related to the neutralino 



sector like the electron electric dipole moment is also interesting 20 
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Appendix 



The contribution to the effective couplings Ql and Qr coming from Figs.(b) and (d) 
are summarized as, 



em, 
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M 4 s/2m w cos (3 



A 



M 4 

2 

+ £^*4 



zf(2,g e jzf(2,g e j + 



2m 2 w cos 2 P 4 \ 



m c2 , A 2 ,„ m„m. 



fie 



(27) 



, , M 2 7 M 4 2m^ cos 2 /5 

where we used the fact m M 3> m e and the kinematical functions F 3 (r) and i*4(r) are 
defined by 

1 



F 3 (r) 
FJr) 



12(1 - r) 5 
1 



-1 + 9r + 9r 2 - 17r 3 + 6r 2 (r + 3) lnr 



1 + 9r — 9r 2 — r 3 + 6r(r + 1) lnr 



(2f 



6(1 - r) 5 

In the limit of — ► 0, -B(/i — > cy) calculated from the photino contribution in Ql can 



be easily checked to be reduced to the result given in ref.]14| 
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Figure Captions 



Fig. 1 

One-loop diagram contributing to the effective coupling Q L of /il £r1- Figs, (a) and 
(b) represent the neutralino contribution and Figs.(c) and (d) represent the chargino 
contribution. There are similar diagrams in which the chirality of the external lines are 
exchanged. Flavor mixings are induced by the off-diagonal elements of slepton mass 
matrices which are expressed by •. It should be noted that the chirality flip occurs on 
the internal line in (a) and (c) and on the external line in (b) and (d). 

Fig. 2 

B(/jl — > e^f) as a function of u in the 77-model. The vertical axis Br stands for 10 11 x 
B(/i — > cy) and the horizontal axis u' should be understood as u = 50 (V + 2). Each line 
corresponds to the various parameter settings for (Mw, sinx) and their values are taken 
as A(80, 0), 5(180, 0), C(80, 0.3) and 5(180, 0.3). 

Fig. 3 

The ratio it! of B(/j, — > cy) against the MSSM as a function of u in rj- model. R is defined 
as R = B V (M W , sin x)/B MSSM (M w ). Each line corresponds to the various parameter 
settings for (M w , sinx) and their values are taken as A(80, 0), 5(180, 0), C(80, 0.3) 
and 5(180, 0.3). 
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